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fT^ ' The de Haas - van Alphen effect in two-dimensional (2D) metals is investigated at 

different conditions and with different shapes of Landau levels (LLs). The analytical 
calculations can be done when many LLs are occupied. We consider the cases of fixed 
particle number (A'^ = const), fixed chemical potential {fj, = const) and the intermediate 

flj ■ situation of finite electron reservoir. The last case takes place in organic metals due to 

r^ ' quasi-one-dimensional sheets of Fermi surface. 

" We obtained the envelopes of magnetization oscillations in all these cases in the limit of 

low temperature and Dingle temperature, where the oscillations can not be approximated 
by few first terms in the harmonic expansion. The results are compared and shown to be 

^ , substantially different for different shapes of LLs. 

The simple relation between the shape of LLs and the wave form of magnetization 
oscillations is found. It allows to obtain the density of states distribution at arbitrary 
magnetic field and spin-splitting using the measurement of the magnetization curve. 

The analytical formula for the magnetization ai n — const and the Lorentzian shape 
of LLs at arbitrary temperature, Dingle temperature and spin splitting is obtained and 

Q ■ used to examine the possibility of the diamagnetic phase transition in 2D metals. 
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O ! I. INTRODUCTION 

f^ . The quantum oscillations of magnetization and conductivity in magnetic field (the de Haas - van Alphen 

f^ ' and Shubnikov - de Haas effects) give abundant information about electronic structure of metals 01 . The 

^^ . theory of these effects in three-dimensional (3D) metals has been constructed g] and proved well by 

many experiments ||l|. However, this theory is not applicable in two-dimensional (2D) compounds, when 

warping of the Fermi surface(FS) is less than the Landau level separation hujc- The calculation of the 

Ch ' kinetic properties of 2D compounds in magnetic field is difficult even in the simplest model approximations 

r^ H] I Hi • Including of the electron-electron interaction, other parts of FS or different scattering mechanisms 

P^ ' makes the problem formidable and the quantitative theory of the 2D Shubnikov - de Haas effect does not 

O > exist. The theory of the do Haas - van Alphen (dHvA) effect is much simpler because the magnetization 

J-^ ' is a thermodynamic quantity. The complete description of the 2D dHvA effect, however, does not exist 

either. There are two approaches to this problem. The first is to make the harmonic expansion of 

magnetization oscillations ^, in analogy to the Lifshitz-Kosevich formula Q. This expansion assumes 

the chemical potential to be fixed (/x = const). In most compounds it is the electron density which 

C^ ' is constant and the chemical potential is an oscillating function of the magnetic field. Even at fixed 

chemical potential this approach in pure metals at low temperature gives the infinite harmonic series for 

magnetization oscillations with slow convergence which is not convenient for the practical use. This series 

has been summed explicitly only at zero temperature and the Lorentzian shape of LLs Iq] . 

Another approach is the direct summation over Landau levels (LLs). It is applicable to the systems 
with fixed electron density (N — const) and just at low temperature T <^ TilJc {TilJc — fieB /m*c is the 
LL separation) since then only few LLs near FS are important H. The analytical calculations of the 
magnetization and its envelope by this method at low but nonzero temperature have been done M, Isl 
only for sharp Landau levels (LLs). This approximation has a very narrow domain of applicability TiiOc ^ 
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T ^ y/ThcUc because the Landau level broadening F has a very strong influence on the magnetization 
oscillations. Any electron reservoirs are also ignored in these calculations although in organic metal 
this reservoir of field independent electron states exists due to quasi-one-dimensional parts of FS. The 
influence of this electron reservoir on the magnetization oscillations has been studied recently |Q| only 
for sharp LLs. Numerical calculations are possible jlQ], O], but the numerical analysis is not evident 
and convenient for practical use. So, in many cases one continues to apply the Shoenberg's 2D harmonic 
expansion. 

Our work is intended to study the 2D dHvA effect with arbitrary shape and width of the LLs and 
under conditions close to real experiments. 

In sections II, III, IV we consider the case of fixed electron density using the method of direct summation 
over LLs. The reservoir electron states and finite temperature are also taken into account. In Sec. II 
we derive the equation for the chemical potential and the relation between the magnetization and the 
chemical potential that are valid at arbitrary shape and width of LLs. In section III this is used to 
derive the simple relation between magnetization oscillations and the density of states (DoS) function. 
We describe in details how this relationship enables the direct measurement of the DoS distribution. 

The problem of the density of states distribution in 2D metals in magnetic field is an important and 
old one. Since the discovery of the integer quantum Hall effect many theoretical approaches have been 
developed to solve it (for example, Q, Q) and it would be very interesting to compare the results of 
these calculations with the experimentally measured DoS distribution. The DoS function tells about 
the electron spectrum and the scattering mechanisms that determine the kinetic phenomena but it is 
difficult to explicitly relate the DoS distribution to the conductivity. The magnetization gives a quite 
direct information about the DoS and an accurate method to extract the DoS distribution from the 
magnetization oscillations is needed. The harmonic expansion is not good here for the following reasons. 
To obtain the DoS distribution function one should consider the case of constant electron density and 
low temperature when many harmonics are important. In this case the explicit harmonic expansion 
(with many terms) of magnetization oscillations does not exist and the influence of LL broadening on 
the harmonic damping can not be separated accurately from other factors. The method of extracting the 



DoS distribution from the dHvA measurements has been proposed by M. Ya. Azbcl |12|. This method 
is valid only in the extremely strong magnetic field when the magnetic length Iq — {cfi/eB)^''^ is much 
less than the impurity Bohr radius. His procedure assumes that one can realize the situation when only 
one LL is occupied (this means that htOc/'^ > Ep , Ep is the Fermi energy). Such a condition can not 
be fulfilled in organic metals and many other compounds where Fermi energy is quite great. Azbel's 
procedure also assumes zero temperature, no electron reservoir and no warping of Fermi surface. Our 
approach has no these defaults and limitations. 

In Sec. IV we obtain the analytical formulas for the envelopes of the magnetization oscillations. In Sec. 
V we study the case of constant chemical potential fi — const using the harmonic expansion. We obtain 
the analytical formula for the magnetization with Lorentzian shape of LLs for arbitrary temperature. 
Dingle temperature and spin splitting. This formula is used to examine the possibility of the diamagnetic 
phase transition in 2D metals. We also derive the relations between the oscillations of the magnetization 
and the oscillations of the chemical potential (in the case N — const) or the particle density (when 
ji = const) at general conditions and find their domain of applicability. 

Both the conditions N = const and fi = const take place in real compounds that was determined on 
the experiment using the shape of magnetization oscillations [|3| , [|l^ , ||lj| or by the observation of /3 — a 
frequency [|5|, jl^l in the presence of the magnetic breakdown. Therefore we consider both these cases. 

II. DERIVATION OF MAGNETIZATION IN THE CASE N=CONST WITH FINITE 

ELECTRON RESERVOIR 

The condition of the constant particle density (N = const) means that the sample is isolated from 
electron reservoirs and the sample volume is constant. This takes place in heterostructures [Q and is 
usually more realistic than restriction of the fixed chemical potential (/i = const). Nevertheless, in many 
solids the exact conditions N = const or fi = const are not fulfilled. The intermediate situation takes 
place, for example, in synthetic metals like intercalated graphite compounds where the magnetic field 
dependent charge transfer between the graphite layers and intercalated layers takes place, and in organic 
metals due to the additional quasi-lD parts of FS, that does not contribute to the magnetization directly 



but serve as a finite electron reservoir. The influence of this electron reservoir on the magnetization 
oscillations has been studied recently |Q] only for sharp LLs and in the main part numerically. We shall 
consider arbitrary shape of LLs. 

Since the total particle number remains constant, this intermediate situation can be considered as the 
case N = const with some additional density of states (DoS). This reservoir part of the DoS does not 
oscillate with magnetic field and changes on the scale of Fermi energy. Since only several LLs near FS 
are important for the dHvA oscillations and the number of occupied LLs in organic metals np 3> 1 this 
reservoir part of DoS may be considered as constant. For the same reason we assume the LL degeneracy 
be equal for all LLs. Then the total density of electron states may be written in the form 

oo _. 

p{e,B)=g{B)Y,Dn{e-nu;,{n+-)) , (1) 

n=0 

where the LL degeneracy g{B) = 2_B/$o + "i? ■ ^c includes two spin polarizations and the reservoir 
density of states nn — const, $o = ^irfic/e is the flux quantum. Dn{s) is normalized to unity and gives 
the shape of the n-th LL with some spin-splitting and the reservoir part of the DoS. 

To calculate the magnetization we need first to calculate the chemical potential. It is given by the 
equation 

r p{e,B)ds 
J l + exp(^) ^ ^ 

It is convenient to count the energy and the chemical potential from the last occupied LL: 

6n = n~ huJcinp + 1/2) ; E = e~ hcjcinp + 1/2) (3) 

where np =hit[Ep/fiujf,] is the number of completely filled LLs when 6fi = . It is an integer and jumps 
by unity one time on each dHvA period. 

Now we apply the direct summation over the LLs. Another approach of the harmonic expansion of 
the DoS (ll|) will be used in the next section for the case ji = const. At T <C hcvc only three LLs near FS 
produce the temperature dependent contribution to the oscillations because the contribution from other 
LLs is small by a factor exp^{—huJc/T). The equation (||) then simplifies 

r3hu>j2 D{E) dE 

n+l= / ^^ ^ (4) 

J-3hujj2 1 + exp ( -^^ 



where h = N/g — np — F/B — np is the filling factor of the last occupied LL, N is the area electron 
density, F = const is the dHvA frequency. n{B) is an oscillating function of the magnetic field. On 
each dHvA period it monotonically decreases from 1 to with increasing of the magnetic field. The DoS 
function D{E) here is periodic D{E + huc) = D{E) and normalized to unity 

hi^c/2 

D{E)dE = 1 (5) 

^hLj^/2 

The LLs just above and just below the last occupied LL arc important when the chemical potential is 
situated just between two LLs {\6fi\ — huc/'^ ~ T). Then these two LLs are in equal position and hence 
have the same shape. Therefore, we have written the same DoS function D{E) for all three LLs. 
Now we have to calculate the thermodynamic potential. By definition 

n = -fcBT /"dep(e,B)lnM+exp^-=^j , (6) 

where the density of states is given by (0). As in the derivation of (H) we calculate the temperature 
dependent contribution only from three LLs near FS. Actually we could leave only one LL but now the 
derivation is more general. The expression for the thermodynamic potential now becomes 

"-^^ H - Tiujcjn + 1/2) 



dED{E)\n(l + e^vJ^-^] 

^3huj2 \ J / 



l''ihuJc/2 

The summation in the first term is easy and gives ghujcinp — 1)^/2 — gf^inp — !)• We need the free energy 



F ^ n + fiN ^ gfi{n + I) + ghuj^ 
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-Tgj^^^ dED{E) In (^1 + exp %^j (7) 

To obtain the magnetization M = ~dF/dB one should differentiate all the magnetic field dependent 
quantities in this expression. As a result one obtains a quite huge formula. We do not write down it here 
because it is useless for the analytical calculations. Even for numerical calculations it is simpler to use 
([7|). The elegant expression for the magnetization can be obtained if many LLs are occupied {np ^ 1). 
This condition is always satisfied in organic metals where np ~ 300. Then in the third term of (0) one 
can differentiate only oscillating quantities since they give additional factor np . The only such a quantity 
is (5/z. We assume the function D{E) to be independent of the position of the chemical potential, that 
is usually well satisfied. The terms with -j^ in dF/dB cancel each other after using (||) and (^. As a 
result we get 



M kC 



1 . Sn 

n + — ^ 

2 hiVr 



(8) 



where S^ is the solution of (Eh. The prefactor 



C = —-hurTip « ^—Ep = const (9) 

B B ^ ' 

is the same for different dHvA periods since g ^ B. The first two terms in square brackets in dq) give the 
sawtooth form of magnetization, and the last term j-^ determines the damping of the oscillations due 
to finite temperature and the LL broadening. The expression in square brackets of (ph is the oscillating 
part of the chemical potential divided by hujc- This can be easily shown using the definition (0) of (5/i. 
Then the formula (0) becomes 

M = C(^ - Ep)/hujc (10) 

Another (more apparent) derivation of this relation between the oscillations of the magnetization and the 
chemical potential (formula ( pq ) ) and its discussion will be given in section IV. 

So, in this section we have derived the formulas (Q) for the chemical potential and (H) for magnetization 
at finite temperature and arbitrary DoS distribution. The density of electron reservoir states was assumed 
to be approximately constant on the scale of huJc- In the next two sections we shall use these formulas to 
obtain more concrete results that are useful for the analysis of the dHvA effect. 



III. RELATION BETWEEN THE MAGNETIZATION AND THE DOS FUNCTION 

Now we derive the relation between the wave form of the the magnetization M{B) and the DoS 
distribution D{E), that allows the direct measurement of the function D{E) at different magnetic field 
values and spin-splitting if the temperature smearing is much less than the LL broadening. We assume 
T ^ r, where F is the width of LLs. Since we need the function M{B) on one particular dHvA period 
(0 < n < 1) it is convenient to consider magnetization as a function of n « np(Bo — B)/Bo, where on 
this period n{Bo) = 0. Differentiating (||) we obtain 



dM _ 
dn 



1 dS^x 
fkOr dh 



(11) 



Differentiating the equation (Eh with respect to (5/i we get 



dii 



3huic/2 



D{E) dE 



?,nu:j2 ^T cosh^ (^^\ 



(12) 



At T <C r the important region of the integration is \E — 5ijl\ ~ T and the function D{E) can be expanded 
in Tailor series near this point 



dSfj. J_ shu,^ 



4T cosh' 



(^) 



-dE 



(13) 



The all odd terms drop out because the integrand should be a symmetric function of {E — (5/i). We shall 
keep only the T^ term. The integration can be extended to infinity. After the integration we get 



^^D{S^,)+T'D"{S^,)^+0{T^) 



Substituting this into (O) we obtain 



dM 
dh 



= C 



-1 



hujc {D{5fi) + T^D"{6^i)^ + 0{T*)) 



(14) 



(15) 



The temperature dependent terms in (nS) are small and can be separated by taking the measurements 
at several low temperatures and extrapolating to T = 0. So, one can measure the function 



dMr^oin) 
dh 



C 



-1 



hujcE>{S^) 



(16) 



On the experiment it is very difhcult to obtain the proportionality coefhcient between the measured 
signal and the magnetization. For this reason we shall look at the constant C as on an unknown factor 
and later we shall describe how one can obtain it. Now suppose one can measure the ratio 
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With this function one can rewrite the equation (fid^ as 

DiSfiih)) 



1 



^-c(^ 



(17) 



(18) 



The function D{5iJL{n)) is not D{5ii) — D{E). The dependence 5^{n) is simply related to the function 
R{fi) since the formula (ll^) can be casted into 



5^{n) = hojc 



R{fi) 



(19) 



The DoS distribution D{E) is just the plot of D{n) (|8|) as a function of (5^(n) (|ig|). 

Now we have to say how one can get the constant C by which the measured torque should be divided to 
obtain the function R{h) (O). The normalization condition (ra) does not give this constant since for any 
C the described procedure gives automatically normalized DoS distribution. The constant C determines 
the strength of the magnetization and, hence, the DoS oscillations. The smaller one takes the constant 
C in ( |l7| ) the larger oscillations of the DoS function he obtains. There is some critical value Cq such that 
if one assumes C < Cq he gets a singularity in the DoS, calculated from dlq). At C = Co this singularity 
is in the middle of a LL. The peak value D{Q) of the DoS is indeed large compared to the average DoS 
D = 1/hujc for small LL broadening T <C fiujc- The correct value of the constant C is then close to Cq 
which gives the function D{E) accurately always except the vicinity of the peaks of LLs. We shall use 
this singularity to obtain more accurate value of the constant C. 

The derivative (Wm has a sequence of the periodically situated minima that occur on each dHvA period 
when the chemical potential crosses the position in the middle of a LL {Sfi = 0). Since np ^ 1 and the 



magnetization is measured on many dHvA periods these minima form a smooth function of magnetic 
field 



MLuiB) = C 



I ' 



hLUcD{0) 



This function monotonically decreases to some finite hmit since huJcD(0) increases with the increasing 
of the magnetic field. This is because fiujc ^ B while the peak value of the DoS I?(0) has the very 
slow dependence on the magnetic field. In the vicinity of the peaks the DoS distribution may be given 
by only two parameters: the width of the LL F and the peak value D{Q) — 1/aT. The width F is 
determined mainly by the impurity scattering which is approximately independent of the magnetic field. 
The parameter a '^ 1 depends on the DoS distribution but also is almost independent of the magnetic 
field. For Lorentzian shape of LLs a = tt. 

Hence, one can assume -D(O) = I/aF — const and from the curve of the minimum values of 



NC„,{B) = C 



flLOr 



(20) 



one can easily obtain two unknown constants: aF and C . This value of the constant C is accurate enough 
to calculate D{E) even in the vicinity of the maxima. 

The above procedure gives the DoS at the position of the chemical potential while the magnetic field 
and the chemical potential are varied. In general, this is not the same as if one fix the magnetic field 
and the chemical potential but vary the energy itself. In our analysis we disregard the dependence of 
D{E) on the position of the chemical potential and hence this difference. This approximation works 
when the DoS distribution is determined by one-electron processes. These are electron scattering on 
lattice imperfections and inhomogeneities, the finite probability of the interlayer jumping and so on. 
When the many particle effects play an important role (for example, the change of the magnetic field 
drives a sequence of phase transitions as in fractional quantum Hall effect), our results are not applicable. 
We have neglected the nontrivial many body effects at the beginning when have written the formulas 
( 2 ) and (^ assuming that the system may be described by the distribution of the single fermion states 
(L). Nevertheless, this has a very wide application region because usually fine many particle effects are 
damped by impurities, finite kz dispersion and other factors, especially when many LLs are occupied. 

One can not get the value of LL separation fiojc without performing some more measurements because 
after we have neglected the temperature dependent terms, the LL separation becomes the only parameter 
with the dimensionality of energy. To find it one needs to compare it with some other energy parameter. 
If one knows the g- factor of the spin-splitting in this compound he can compare UliJc with the spin-splitting 
energy. One can also use the temperature dependence of the amplitude of magnetization oscillations to 
obtain fiuJc , but this would be less accurate than the cyclotron resonance method. So, there are many 
ways to get TujJc- 

To summarize, the suggested procedure of extracting of the DoS distribution from the oscillations 
of magnetization consists of the two steps. The first is to obtain the constant C (that normalizes the 
measured signal) and the second to plot D{E) itself. To do this one should measure the torque as a 
function of x = F/ B where F = const and x changes by unity on each dHvA period. Then x — np -f n- 
This function M{x) should be measured at several low temperatures and extrapolated to T = 0. As 
a result one gets an oscillating function Mq(x) with period equal to unity. We also need its derivative 
Mq{x) — dMQ{x) / dx which is also an oscillating function. The minima of AfQ(a;) form a smooth monotonic 
(x). Extrapolating this function to a; = one gets the constant 

C = hm Mjnin{x) 



Or one can use the formula (20) and obtain the constant C from Mmin{x) more accurately. Substituting 
the functions R{x) — Mq{x)/C and R'{x) = Mq{x)/C into (llq) and (noh one can plot the expression jlq ) 
as a function of (na) on one period of the oscillations. This plot is the desired DoS distribution at the 
magnetic field corresponding to the chosen dHvA period. 

The described procedure allows to measure the DoS distribution on the Fermi level for different mag- 
netic field values (because one can obtain D{E) on each dHvA period) and for different spin-splitting 
energy (because one can tilt the magnetic field with respect to the conducting plane of a sample) . This 
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information about the DoS distribution is very important for the study the role of different scattering 
mechanisms in the electron motion at different external parameters. Such an information can not be ob- 
tained using the previous methods of processing the data of the dHvA measurements where one assumed 
some particular shape D{E) and then numerically calculated the magnetization. The results of the cal- 
culation (usually only of the envelope) one compared with the measured signal and if the agreement was 
good enough the chosen function D{E) was taken as a result. Our method is more simple and accurate. 

IV. THE ENVELOPE OF MAGNETIZATION OSCILLATIONS 

We can also calculate the envelope of magnetization oscillations for several simple DoS distributions 
in the limit hujc ^ F 3> T. At low temperature the envelope turns out to depend strongly on the shape 
of the LLs while the additional constant contribution to the DoS from other parts of FS almost does 
not change the envelope. The envelope of magnetization oscillations does not give as much information 
about electronic structure of the compounds as the DoS distribution does. Nevertheless it is still useful 
for the analysis of the dHvA effect and is quite easy for measurement. 

To calculate the envelope of magnetization oscillations we shall take the total density of states function 
in the form 

where F is the width of LLs, k < is a number that determines the constant part of the DoS and 

n^./2 fE\dE 



converges rapidly. Then the function D{E) is also normalized to unity. The degeneracy g{B) of LLs 
should be renormalized to include all additional parts of FS. 

How does the envelope of magnetization depend on k. The answer is that in the limit np 3> 1 the 
constant part of DoS affects only the shape of magnetization but not the envelope of the oscillations. This 
is different from the case /i — const where the constant part of DoS does not change both the envelope 
and the shape of the oscillations. 

To show this we first should substitute ( [2l| ) into (H). The LLs just above and just below the last occupied 
one contribute only when \5^\ « TiliJc/2. As we shall see, if T, F <C fioJc the extrema of magnetization 
take place when \5^\ ^ ?iWc/2. So, the regions where {fiujc/'^ — |<5/i|) ~ T ^ hujc/2 are not important for 
the envelope of magnetization and T and F dependent contribution is given by only one LL. Then the 
equation (0) becomes 



hLJc/2 



D{E) dE 



n<^o/2 i + exp(:^) 



(22) 



The integral 



Sfi 



and after substitution of (^21) the equation (^23) acquires the form 



huj^/2 



{1- k)Do{E) dE (\ Sfi 



K - 



Substituting this into the expression (|^) for the magnetization we get 

'1 Sn 



(23) 



9* 

M = -—haJrUf 

B 



2 ' huJr 



(24) 



where g* = {1 — K)g is the LL degeneracy of only the oscillating part of DoS and 

M^N ^ / /^ . N (25) 

The magnetization on each dHvA period may be considered as a function of (5/i : M{B) = M{6fi). The 
envelope of the magnetization oscillations is then given by 



M±{B) = M{5ii,^{B)) = C* 



-X + -T no{d^iex) 

A nujc 



(26) 



where C* — g*huJcnF/B and the extremum values Siiex{B) at which the magnetization has maxima or 
minima are given by the equation 

dB dSfi \1ii^c dSjiJ 

After accounting for ( p3|) this equation becomes 

dh _ [^'^'''^ (1 - K)Da{E) dE 



dSn J-hu.j2 ATcosh^ ( '^-^pA 



huJc/'^ 



Do{E) dE 



-fic^e/2 4TcoshM^^^ 



K 


1 


huJc 


TlLJ 


1 
TthJc 





(28) 



The equation (28) for S^^x is independent of k and so does the envelope of magnetization given by (Eq). 
The function Sfi{B) is monotonic on each dHvA period and different on different dHvA periods. More- 
over, if D(E) is a symmetric function then Sfj, is an antisymmetric function of {n — 1/2) that can be 
obtained from (0) using the identity 



1 _ 1 1 sinh f - sinh ^ 

1 + exp -^ ^ 2^2 cosh f + cosh ^ 



The magnetization (H) then is also an antisymmetric function of {h — 1/2). 

In the limit T ^ F the function Dq{E) may be expanded in Tailor series as has been done in (p^). If 
we rest only ^ T^ terms the equation (pSl) becomes 



7^2 1 

Do(^Mex) + T''D'^{5^l,,)— = -— (29) 

6 TiLOc 

Formula(^5[) can also be simplified. Integrating by parts and expanding up to the terms ~ T^ we get 

no = Go{5^i) + T^'D'^{5^Ji)'^ (30) 

where we have introduced the function 

/E 
Do{E')dE' (31) 

-huJa/2 



The function Go{E) is dimensionless and changes in the range [0;1]. From the equation (|28| ) one can 
immediately say that if Dq{x) has the exponentially falling tails, the 5^ex ~ ±r hi{'hu)cjT). If Dq{x) falls 
off as x~^, then 5^ex ~ F(?ia;c/F)^/'''. This determines the behavior of the envelope (pq). Now we shall 
consider in details two different functions Dq(x). 



A. Exponentially falling D(x) 

There are many symmetric exponentially falling functions that are eligible to be a density of states 
function Dq{E). They all lead to the close expressions for the envelope of the magnetization oscillations 
and for example we shall take 



1 

4rcosl?(^) 



D^iE) = .^^^,,,^, (32) 



Now the equation (Eq) possesses the symmetry T ^^ F if T or F <C fiuJc and one can easily obtain the 
limit F < r from the limit T < F. 

The equation (E3) can be solved by the iteration procedure with small parameter T^/F^. In zeroth 
approximation 

^o(<5a.L) = Ijj-r^ = ^ 

4Fcosh^[%^j ^'^'= 

Since the magnetization m) is an antisymmetric function of 5fi we shall consider only one (negative) root 
of this equation 



<5^L = -2rarccoshJ^ « -Fin (^ ) (33) 



From equation (au) we get 



"o(^M°J « — ^ (34) 



Substituting (B3h and (M) into ( |26| ) we get the envelope in zeroth approximation 

2 hcoc V F / nujc 



(35) 



If one makes the replacement T —^ T this result coincides with the result of Vagner et al. Q , obtained 
for the case F = and finite temperature. It's not surprising because of the mentioned symmetry F <-> T 
with the DoS function (p^). 

One can easily obtain the first temperature correction to the envelope (^5|) . In the first approximation 
the equation (|2^) becomes 

DoiSfie.) + T'^D'iiSy^l,)^ = ^ (36) 

The correction to the chemical potential 

A — ^ 1 JC ^ ^Q\"rex) 6 /o'7\ 

A^Me. = ^/.e. - Sf,,^ = ^,(^^oj (37) 

This correction Aj^/i ~ T^. The first correction to the magnetization AA/+ ^ (Aj^/i)^ because according 
to the equation G^ M' {S^^^) = and the correction to the envelope of magnetization 

d^M{S^l%J (ATA^e.)^ _ . d^no{S^^°.) (ATMe.)^ , . 

rf('5MOJ2 2 d(<5/.oj2 2 ^ > 

In writing this we have used M'{5^1^) = Af'((5/ie^)— M"{6fj,'^^)ATfiex- From (pO|) we have 



n'i{S^^l) ^ D',{S^^'i, 



(39) 



In our case (P2 



D'oiSfil) = - 



sinh 



2r 



4r2 cosh 



\ 2r 



ThuJr 



msf^'^j 



3 . ^ T'^' 

2f5^ and A^^e. = --- 



Substituting this into (|3q) and (B9|) and using (p5|) we obtain the envelope in the first approximation 







1 


r 




^n^c^ 


r 


T 


1 


^7r2^ 


2 


/J^\ 


4" 


± 


= c* 






In 





















2 


hujc 




V i' J 


TlLJc 


hujc 


2 


V4, 


/ 


U'y 





(40) 



B. Lorentzian shape of LLs 



For the Lorentzian shape of LLs 



Do{E) 



I/ttF 



1 + {E/ry 



(41) 



we shall do the same steps as for the exponentially falling D{E). From equation (E9[) in zeroth approxi- 
mation we obtain 



hUJr 



5^il,^T\|-^-l^.Jvn^ 



and 



no{5^ll^ 



■ arctan ' '^''^ 



Substituting this into the expression for magnetization (pq) we get 



Ml = ±C* 



1-2 

2 V irhidr 



(42) 



(43) 



(44) 



This is different from (pq). So, the envelope is different for different shapes of LLs. The temperature 
correction is also different. 

To obtain the first temperature correction we need 



D',{5nl 



'"'■D 



hujc V Thujc 



and D[^{Sn[ 



ex) 



dfK 



T{huocf 



^TtJ'e 



7r5/2T2 



2^/TfwJc 
Substituting this into (|3§) and (^ and using (Q) we get 

M^ = C* "^ 



2 V TThuJr 



^brp4 



ttHuJc 8T'^{fiuJc)'^ 



(45) 



The formulas (HO) and (Hq) are not very confident even at low temperature. This is because the 
oscillating part of the DoS distribution function depends not only on the one parameter - the width F, 
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but also on the LL separation hujc and hence on the magnetic field. Any case the approximate functions 
( p2| ) and (^l]) may be too far from reality. The exponentially falling Dq{E) is always the incorrect 
approximation, while the confidence of the Lorentzian shape of LLs in 2D metals will be discussed later. 
The envelope of magnetization oscillations is very sensitive to the tails of DoS function D{E) at low 
temperature and one should use more realistic (probably magnetic field dependent) functions D{E) for 
accurate calculations. 



Nevertheless, the formulas (40) and ( }45| ) can be used to indicate the qualitative features. They predict 
that the first temperature correction to the envelope is proportional to T^ and that this envelope depends 
strongly on the shape of LLs at low temperature. 

We can now say also at what conditions the previous results are valid. At TV = const and finite 
temperature the only results m , B for magnetization and its envelope have been obtained for sharp LLs 
{D{E) — S{E)). These formulas have been derived also in the limit T <C hujc- If we assume the shape of 
LLs to be Lorentzian these result are valid if huJc ^ T ^ ^/Tfujc- This narrow region may not exist at 
all, that explains why many experimental data can not be described by these formulas. 



V. MAGNETIZATION OSCILLATIONS AT CONSTANT CHEMICAL POTENTIAL 

The case of constant chemical potential may take place in organic metals when the electron reservoir 
formed by the one-dimensional parts of FS is very large. This has been observed recently on /3"-(BEDT- 



TTF)2SF5CH2CF2S03 |14| by the inverse sawtooth wave form of the magnetization oscillations. 

The case /i = const may take place also due to giant magnetostriction (volume change) that cancels 
the variations of chemical potential associated with the oscillations of the density of states. This has 
been first observed in beryllium [Q . In organic metals the electron density is not usually large enough to 
observe the giant magnetostriction. The restriction /i = const can be achieved, in principle, by connecting 
the sample with some big 3D piece of metal. 

The case fi = const is interesting because it is more favorable to observe the diamagnetic phase 
transition in 2D metal and Condon domains pi] ]. It is also interesting to compare the cases /i = const 
and N = const. 

Although in pure 2D materials at low temperature the dHvA oscillations are essentially non sinusoidal, 
the expansion over harmonics is still useful and the infinite harmonic series for some special densities of 
states can be calculated with good accuracy. The 2D harmonic expansion of the magnetization oscillations 
has been first obtained by Shoenberg ||5| in the case fi = const using the intuitive phase smearing 
arguments. In the special case of the Lorentzian shape of LLs it has been derived recently ||^ using 
Green's function. We shall give another simple derivation of this harmonic expansion for arbitrary shape 
of LLs. By definition the thermodynamic potential is given by the same expressions (^,|l|) as in the case 
N = const. We assume D{x) to be independent of the number of a LL, since in dHvA effect only the 
LLs near FS are important and for these levels the function D{x) does not differ substantially. Then 
p{E) = p{E/fiujc) is a periodic function. One can take the spin-splitting into account by introducing 

Ps{E) = p{f^ ^/2) + p(|[5~ '^ ^/2)i where A is the energy of spin splitting divided by huJc. This is 

equivalent to introducing fls{p) — ^^(7-^ — ^/2)+^(7^ + A/2), since integration in (^ is between infinite 

limits. Spin wiU be recovered in the final result by introducing Ms{p) = M{j^ A/2) + M{jJfj- + A/2) 

, where M{fj,) is the magnetization without spin. 



The density of states (pj) can be rewritten as a harmonic series 



k— — oo 



2TTt ( \ k 

riuor 2 



Mk)l±I^ (46) 



where we introduced the Fourier transform of the DoS distribution on each LL 

A{k)^ I ' D{E')e-x^(2m$^dE' (47) 

J^n^ \ huJcJ 

and the factor (1 + sign(i?))/2 indicates that no electron states exist for _E < 0. One can omit the term 
with fc = since it doesn't affect the magnetic oscillations. Substituting (Eq) into m) and integrating 
two times by parts we get 
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27ri ( i 1 fc 

TiLUr 2 



A{k) (48) 



Since §^ ^ only in a small region near FS, the factor (1 + sign(_E))/2 has been omitted. 
Introducing the Furrier transform of -g^ 

FT(k) .JdE ^-^j exp (2..L_^ j = ^,^,(,^,4^,^^) - -^ (49) 

we get 

Differentiating (pO[) we obtain the harmonic expansion of the magnetization oscillations 



M = -g |„.= fl f tip 3in (£^2..) m FAk) (51) 

fc=l \ c / 

This coincides with the Shoenberg's formula Q. 
The oscillations of the electron density are given by 



TV = -^ \B=const= J J2 ^T~ '''' iiir^'''') ^^^^ ^^^^^ ^^^^ 



9r2 

9^ 



fe=l 



From (pnj and (p3) the simple relation follows 

M = -iiN/B. (53) 



It is valid only in the limit /i 3> fiiOc where the formula ( |5l| ) is valid. 

In the case N = const one can obtain the similar relation between magnetization and the oscillations 
of the chemical potential. The magnetization oscillations at A^ = const 



M 



d{n + Nn) _ dn 



N— const — rs j-y Ifi^N— const 



dB I"—"- QB 



_ dn at] ^\ -_^ 

I o \N ^B— const 1-^^ I J j-^ \N— const no \fi^N— const 

So, the formula ( pl| ) is valid also in the case N — const but the chemical potential becomes an oscillating 
function of the magnetic field. The total electron density 

N^ [noiE) + nR{E)]dE + N ^ const (54) 

-' — oo 

where no{E) — {g/hu!c)0{E) is the average DoS on the 2D part of FS and nu{E) is the reservoir DoS. 
Substituting N = J^^no{E) + nR{E)]dE and (|5|) we get 

M = ^ r [n,,{E) + nR{E)]dE (55) 

^ JEp 



If nR{E) ~ const this simplifies 



M^^{no+nR){^i-EF). (56) 



This expression coincides with the formula (H). The variations of the chemical potential \ii — Ep] < ?iWc/2 
and the relation ( |5^ ) is valid if n'j^{EF)'huJcl'2. <C nRlEp)- This is usually fulfilled in the limit ^/fiuic '> 1. 
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The formula (|5^) has been obtained previously only for sharp LLs in |9|] and for Lorentzian shape of 
LLs without reservoir states in l6l . We have derived this formula in general case so that its applicability 
region becomes clear. 

The electron reservoir (the independent of the magnetic field DoS) does not contribute to the mag- 
netization oscillations ai n = const. This is evident from (ph since fl is linear in p{E, B) and the DoS 
p{E,B) is the only magnetic field dependent factor in right hand of the expression (0). So, in the case 
fjL = const one can consider only oscillating part of the DoS. 

The expression ( [si] ) depends essentially on the function A{k) that is determined by shape of LLs. The 
question about the DoS distribution on the LLs in 2D metals is open. The Lorentzian shape of LLs (E^) 
at small LL broadening F leads to the Dingle low of harmonic damping 

A(fc).exp(-I|) (57) 

This exponential low of harmonic damping has a direct physical meaning - it assumes the noncoherent 
scattering of the electrons on the impurities. The width of LLs F is related to the electron scattering 
relaxation time r = irU/T and ( pTj ) is the amplitude of the probability that an electron returns to the 
initial state after k cyclotron periods. This can be shown by follows. The DoS is connected with the 
electron Green's function by the relation 

p{E) = - -ImG^ (r, r, E) (58) 

The imaginary part of the retarded Green's function 

1 1 1'°° 

IniG^{r,r,E):^ -{G^ -G^)=-- / i < ip+{r,Q)ij;{r,t) + iP{r,t)^+{r,Q) > e^p{~iEt/K)dt (59) 

^ ^ J -ao 

Substituting this into (^) we get 

A{k) - j p{E) exp f ^^ j dE = 



27r 
using 

we finely obtain 



1 /"^ 1"^ / 2'KiEk iEt\ 
— di < V+(r-,0)V'(r,i) +V(r,i)^+(r,0) > / d^exp — ^]dE 



1 f°° , {2mEk iEt\ , ,, , , , 

— / dE'exp — —\dE^6{t- lixkju:^) 



27r ./_^ V 'hiOr h 



A{k) ^<tP+ (r, 0)'iP{r, T^k) + ^(r, Tcfc)V'+ (r, 0) > (60) 

where Tc = 2t: /ojc is the cyclotron period. 

The formula (pj) states for the arbitrary ratio T /fiujc and is more general than (p|). In 3D metals after 
each scattering process an electrons leave the impurity for a long time because it is free in the direction of 
the magnetic field. In 2D case the electrons can not move in z direction and return to the same impurity 
after each cyclotron period. If the phase smearing time is longer than the cyclotron period the electron 
interfere with itself and this may change the amplitude (p^) and hence the shape of LLs. So, the Dingle 
low of harmonic damping may violate if the cyclotron energy is greater than warping and the width of 
LLs fiLOc ^ VF, F. Nevertheless on many experiments on organic metals no remarkable deviation from 
the Dingle law of harmonic damping has been observed. This may be because even small kz dispersion 
breaks the coherence of 2D scattering. We now assume the Dingle low (|5^) to be valid and try to sum 
the series (|5^) for magnetization oscillations at arbitrary T and F. 

The oscillating part of the magnetization is 
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where 



M-^ImY.eM^yk^m-lr^-^, (61) 

fc=i 



A — 2n'^T/huj , b = tt/uJcT and y — — — (62) 



huj. 
To calculate this sum we expand 

1 2e-''^ 



sinli(fcA) 1 - e-''^ 
Now sum over k is the sum. of geometric series and can be calculated 



= ^2e-'=^e-'P^^ (63) 



p=0 

where 



M=^Imf^Qip) (64) 



nc ^ - exp [iy - b - A(l + 2p)] 

^(^) = l + exp[.y-b-A(l + 2p)] ^'') 



Q{p) is a smooth monotonic function so one can apply the Euler-Maclaurin summation formula p9| | 



oo -, <.oo oo (-)(2j-l)(n\ 



where Bk are the Bernoulli numbers. The terms of the sum in the right-hand part of eq. (pq ) fall off 
~ A. I j^y , hence their sum converge rapidly. The first term ~ ^ Q(0). So, all the sum in the right-hand 
part is much less than Q(0) and we can ignore it. The integral 

/■°°^. ^, 1 f sin(y) exp [-5 - A] \ 
Im / Q(p)dp^— arctan i , , , \, (67) 



Substituting (^5|) , (|6g) and ( |67|) into (|64| ) and recovering the spin splitting we obtain 






sm(y — ttA) exp —6 — A 

arctan | — — ; -r 

1 + cos(y — ttA) exp [—b — \\ 

A sin(y — ttA) exp [—6 — A] 



1 + 2 cos{y - ttA) exp [-6 - A] + exp [-26 - 2A] 



-K{(y-^A)^(y + ^A)} (68) 



This expression gives the amplitude and the shape of magnetization oscillations. It is much more 
convenient for practical use in the limit of weak harmonic damping than the infinite series ( pl|) . In the 
limit of strong harmonic damping A ^ 1 or & ^ 1 formula (pq) correctly gives the first harmonics. In the 



other limit A = the formula (68) becomes exact and simple 



A* / sin(w — 7rA)e ^ \ ,, ., , . ,, 

Ms = -^ arctan ^—, '-— — r + { y - ttA -> w + ttA } 

This coincides with the recent result of T. Champel and V.P. Mineev 0. At 6 = , using the identity 

sin a; \ /sin(x/2) 

arctan = arctan ' 



l-t-cosxy \cos(a;/2) 

one gets a saw-tooth function as has been predicted by Peierls [pO| . 
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If spin splitting energy is just n-times the cyclotron energy (this can be achieved by tilting the magnetic 

(69) 



field) and chemical potential fi ^ hujc , formula (68) substantially simplifies 

2gn r / sinuexp [— 6 — A] \ Asin ycxp [— 6 — A] 

arctan ' ' 



ttB I \1 + cosyexp [— 6 — A]/ 1 + 2cosy exp [— 6 — A] + exp [— 26 — 2A] 

Now one can obtain the envelope of magnetization oscillations. The equation ^- = gives 



-(1 + 3ea;2 + A + A ex^) + Jil + 3ex^ + A + A ex^)^ ~ 86x2(1 + 2A + ex^) , , 

cos y* = ^ '—^ —^ 1 ^ i (70) 

where ex = exp[— 6 — A], that determines the positions of extrema of M{B). Substitution this into (|& 
gives the envelope of magnetization oscillations. In the limit A^l, 5<Cl, X ^ b we have 



M±{B)^±^^ 



2V2& 

$0 



1-^— + 



(Ai,5i) 



(71) 



where A and b are given by (p2) . The result M± does not depend on temperature in first approximation 
at A ~ 6 ^ 1. This is not surprising because A{k) — exp(— fefc) has stronger dependence on b than FT{k) 
(eq. E^) on A for small k, that are important in the sum (|5l|). Since b = 2itV /ftiOc the formula (ItT 



gives the same envelope as (44). Thus, if the Dingle low of harmonic damping is valid, the envelope of 
magnetization oscillations is identical in the cases of constant particle density and fixed chemical potential 
in the limit b — 2-kT /fiuJc ^ 1 and zero temperature. 

Using formula ( pSJ) one can examine the possibility of the diamagnetic phase transition(DPT) Q. It 
must take place if dM/dB > 1/Att. It has been observed in some pure 3D metals due to the formation of 
Condon domains |£1[. In 2D case it can lead to an energy gap at the FS, analogies to that of QHE. The 
diamagnetic phase transition is not possible in available heterostructures since the magnetization is ~5 
orders less than necessary. In organic metals the electron density is much greater than in heterostructures 
and the constraint fi = const is more favorable to observe the DPT. For example, in the compound 
K - {BEDT-TTF)2h at Td = O.IK the slope dM/dB turns out to be only 5 times smaller than needed 
for this phase transition. Maybe there exist another compound where this derivative dM/dB is great 
enough. 

There is another characteristic parameter of magnetization oscillations, besides the frequency and the 
amplitude, that can be compared with the experiment. It is the value of the magnetic field Bq above which 
the spin-splitting is clearly visible on the graph M{B). More formally, it is the point Bq, above which 
the first derivative of magnetization has an additional zero on each period due to the spin-splitting. This 
point Bq is extremely sensitive to the parameters of the system, such as the tilting angle of the magnetic 
field, the shape of LLs and the reservoir DoS, the temperature and so on. But to obtain an accurate 
theoretical description of this critical field one needs to make the numerical calculations to take into 
account many factors, unimportant for the amplitude of the oscillations. 

In the 2D organic metals the magnetic breakdown W usually takes place at some field Bmb- The 
frequency of dHvA oscillations is different below and above Bmb ■ In the vicinity of Bmb one should 
apply more complicated theory of the dHvA effect ( ||ll| and the references therein). But far from Bmb 
all formulas remain valid, but one should always substitute fi by hujc ■ F / B , where F is the new dHvA 
frequency. 



VI. CONCLUSION 

The dHvA effect in normal 2D metals has been investigated at different conditions and at the arbitrary 
shape and width of LLs. Finite temperature, spin-splitting and electron reservoir states are taken into 
account. In the limit Ep/hcuc ^ 1 the formulas greatly simplifies and many results can be obtained 



analytically. In this limit there are simple relations (p3) and (56) between the magnetization oscillations 
and the oscillations of the chemical potential (in the case N — const) or the oscillations of the electron 
density (in the case fi = const). These relations are valid at the arbitrary shape of LLs. 

We have obtained the analytical formulas for the magnetization and its envelope at different conditions 
and shapes of LLs, that may be useful to treat the experimental data. The envelope of the magnetization 
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oscillations depends strongly on the shape of LLs in the limit of weak harmonic damping (compare the 
formulas 1^ and H). 

We derived the explicit relation between the wave form of the magnetization oscillations and the shape 
of LLs at arbitrary reservoir DoS. This relation makes possible to obtain the DoS distribution at arbitrary 
magnetic field and spin-splitting from the measurement of the dHvA effect. The DoS distribution gives 
much more information about the electronic structure of compounds than the Dingle temperature and 
the effective electron mass. This may help to solve the theoretical problem about the DoS distribution 
in 2D metals in magnetic field. 
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